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Shafranov's virial theorem implies that nontrivial magnetohydrodynamical equilib- 
rium configurations must be supported by externally supplied currents. Here we extend 
the virial theorem to field theory, where it relates to Derrick's scaling argument on soliton 
stability. We then employ virial arguments to investigate a realistic field theory model 
of a two-component plasma, and conclude that stable localized solitons can exist in the 
bulk of a finite density plasma. These solitons entail a nontrivial electric field which 
implies that purely magnetohydrodynamical arguments are insufficient for describing 
stable, nontrivial structures within the bulk of a plasma. 
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Ideal single fluid magnetohydrodynamics obeys an integral relation which is known as 
Shafranov's virial theorem |l|, ||. It implies that an ideal magnetohydrodynamical sys- 
tem can not support nontrivial localized structures. Instead any nontrivial equilibrium 
configuration must be maintained by externally supplied currents, a guiding principle in 
the design of contemporary magnetic fusion devices. 

Ideal magnetohydrodynamics is supposedly adequate for describing the ground state 
equilibrium geometry of a plasma. It also provides a starting point for a weak coupling 
Bolzmannian transport theory 0. But as an effective mean field theory it lacks the 
kind of detailed microscopic information which is needed to properly account for the 
electromagnetic interactions between the charged particles within the plasma. For this, 
ideal magnetohydrodynamics should be replaced by an appropriate classical field theory 
of charged particles. With a firm microscopic basis and established set of rules for 
systematic computations, a field theory model can provide a rigorous basis for describing 
thermal fluctuations and dynamical effects, including transport phenomena and issues 
related to plasma stability and confinement. 

Here we extend Shafranov's virial theorem to classical field theory where it yields a 
variant of Derrick's scaling argument, widely employed to inspect soliton stability. We 
then apply virial arguments to a realistic field theory model of plasma. In accordance 
with ideal magnetohydrodynamics, we conclude that the field theory does not support 
localized self-confined plasma configurations in isolation, in an otherwise empty space. 
But when we inspect the finite density bulk properties of the field theoretical plasma, we 
find that the virial theorem does allow for the existence of stable solitons. These solitons 
describe extended collective excitations of charged particles in the otherwise uniform 
finite density environment. Our results are consistent with a recent proposal ||, that 
a finite density field theoretical plasma supports stable knotted solitons ||. Indeed, we 
expect that these solitons can be employed to describe a variety of observed phenomena. 
For example coronal loops that are present in solar photosphere are natural candidates. 
The properties of these solitons may also become attractive in fusion experiments, where 
their stability might help in the design of particularly stable magnetic geometries @. 

Shafranov's virial theorem [1] follows from the properties of the magnetohydrody- 
namical energy-momentum tensor T^ u in the ideal single fluid approximation. Its spatial 
components are 

T* j = pv i v j + (p+-B 2 )5 lj - B l B J (1) 

while the purely temporal component coincides with the internal energy density, 

T 00 = -pv 2 + -B 2 + -?— (2) 

T 2 7-1 y ' 

Here 7 is the ratio of specific heats. The fluid variables are the mass density p, the 
(bulk) fluid velocity v l and the pressure p, and B 1 is the magnetic field in natural 



units with fiQ = 1. The plasma evolves according to the Navier-Stokes equation which 
follows when we equate the divergence of the energy-momentum tensor with external 
dissipative forces. These dissipative forces are present whenever the plasma is in motion, 
but cease when the plasma reaches a stable magnetostatic equilibrium configuration that 
minimizes the internal energy 

E = fd 3 xT 00 (3) 



Shafranov's virial theorem follows when we subject (0) to a scale transformation of 
the spatial coordinates x l — > Ax* with A a constant. For the magnetic field we select 
Bi(x) — > X 2 Bi(Xx), as customary in Maxwell's theory. But for the pressure p ideal 
magnetohydrodynamics does not supply enough information to determine its behaviour 
under a scale transformation. For this we assume that the pressure is subject to the 
standard thermodynamic scaling relation of a thermally isolated gas, 

pV 1 = constant (4) 

This implies that under a scaling p(x) — > A 37 p(Ax). If we assume that the value A = 1 
corresponds to an actual minimum energy configuration of the energy (|3|), when viewed 
as a function of A the energy (|3|) then has an extremum at A = 1. Consequently 



= -nr^ -~ ! dx \ 3v+ 2 B ) s -J dxT '- (5) 

The magnetic contribution to the pressure is manifestly positive definite. Furthermore, 
(collisionless) kinetic theory relates the pressure p to the kinetic energy of the individual 
particles, which is similarly a positive definite quantity. The integrand in fl5|) is then pos- 
itive definite, and we conclude that under the present assumptions non-trivial localized 
equilibrium configurations do not exist in ideal magnetohydrodynamics []J, @. 

We have formulated our derivation of Shafranov's virial theorem so that it relates to 
Derrick's scaling argument in classical field theory |J. For this, we consider a generic 
three-dimensional Hamiltonian field theory model with classical action 

S = f dtd 3 x C(i/j) = f dtd 3 x{7T a ip a - H[ir,(p]} (6) 

The fields ip a ~ (7r a ,ip a ) are canonical conjugates with Poisson bracket 

K(x),yA</)} =SAx-y) (7) 

Notice that the time derivative in ([]) acts asymmetrically. But in the sequel it will be 
useful to consider symmetrized quantities, and for this we generalize the time derivative 
term by a canonical transformation into 

/ dtd 3 x n a ip a — > / dtd 3 x [ aTT a ip a + (1 - a)ir a (p a } (8) 



where a parametrizes the canonical transformation. 

We assume that the Hamiltonian if is a functional of the fields ip = (n, ip) and their 
first derivatives only, with no explicit dependence on the space coordinates x % and time 
t = x . We then obtain the energy-momentum tensor directly from Noether's theorem: 
Since there is no explicit dependence on x M 

!£ = W^ a + i^f a Wa (9) 

and by employing the equations of motion we identify the components of the energy- 
momentum tensor with the ensuing four conserved currents 

T ^ = Jrr d ^" - 5 ^ C ( 10 ) 

In general fllOD fails to be symmetric. But in the following we find it useful to consider 
symmetrized quantities, and for this we can re-define 

T\ -> T\ + d p X"\ (11) 

where X p ^ u = —X w v has no effect on the dynamics. (Note that if the theory fails to 
be Lorentz invariant, in general there will be no symmetry between the momentum flux 
T°i and the energy flux T l .) 

We are interested in a scale transformation of the spatial coordinates x l —>■ Xx l , 
which sends ip a (x) — > \ Da i/j a (\x) . Here D a is the scale dimension of the field ip a . By 
considering an infinitesimal transformation with A = 1 + e we find for the generator 5s 
of the scale transformation 



Ss^a = X l di7l a + 


D>a 


(12) 


5 s ^ a = x l di V a + D^ a 


(13) 


the energy density this yields 




S S T° = ( -T\ + 9 i (x i T° ) + 5> Q 


\5T° Q , 5T% a , " 

XI ^ a + xa i dk ^ a 
_ bip a dd k ip a 


} (14) 



In general the scale dimensions can be arbitrary, and there is no a priori relation between 
the different D a . But if the scale transformation is a canonical transformation it must 
preserve Poisson brackets, which implies that the scale dimensions of a canonical pair 
are subject to 

D a , + D* a = 3 (15) 

The generator 8g of such a canonical scale transformation can be computed from Noether's 
theorem, and by properly selecting the value of a in (§D we arrive at the symmetrized 
form 

§c = fd 3 xx l T°i (16) 



We are interested in a refinement of ideal magnetohydrodynamics, a microscopic 
field theory model of a two-component plasma with negatively charged electrons (e) and 
positively charged ions (i) and classical (first-order) Lagrangian ||, 

c = E k d t A k + l - (r e d t ^ e - d t rje + rAi>i - d t rM - \®i - \b\ 

- 77- \{d k + ieA k )^j e \ 2 - —{(dk-ieAM 2 + A { d k E k - e^ e *^ e + e^*^) (17) 
2m 2M 

Here ip e and ipi are (complex) non-relativistic Hartree-type fields that describe electrons 

and ions with masses m and M and electric charges ±e respectively, together with their 

electromagnetic interactions. Note that we have realized Maxwell's theory canonically 

so that the electric field Ei and spatial gauge field A{ form a canonical pair, with the 

temporal Aq a Lagrange multiplier that enforces Gauss' law. Since the time derivative 

appears linearly in the charged fields, the action (|TT| ) admits a proper Hamiltonian 

interpretation with if)* i the canonical conjugates of tp e ^. Notice that for definiteness 

we have chosen both charged fields to be commuting. This should be adequate in the 

Bolzmannian limit, relevant in conventional plasma scenarios where the temperature is 

sufficiently high so that bound states (hydrogen atom) are prevented but not high enough 

for relativistic corrections to become important. Notice that we have also introduced an 

appropriate symmetrization of the form (|8]) in the time derivative terms of the charged 

fields. Finally, besides the terms that we have displayed in ( |T7D we implicitely assume 

the presence of chemical potential terms that ensure overall charge neutrality. However, 

fr the present purposes such terms are redundant and will either remain implicit, or will 

be enforced by appropriate boundary conditions. 

We propose that the advantage of (|TT| ) over ideal magnetohydrodynamics is, that 
(|17j) provides a firm microscopic basis for systematically computing various properties 
of a plasma. For example an appropriate version of the equation of state (|]) can be 
derived from (|l7D. In particular, (|T7|) yields immediately the standard electromagnetic 
many-body Schrodinger equation for a gas of electrons and ions. 

The energy-momentum tensor T^ v can be computed directly from (p~0|) . After we 
introduce an appropriate symmetrization which ensures manifest gauge invariance, we 
find for the energy density 

T° = ^{sin 2 a|A^ e | 2 + cos 2 a\D* k ^\ 2 } + — + — - A (8^ + e[<^ - VCVeD 

(18) 
where D k = d k + ieA k and /i = m sin 2 a = M cos 2 a is the reduced mass. For the 
spatial components of the energy-momentum tensor we find similarly, with the help of 
the equations of motion 

T\ = EiE k + B z B k - ^jsin 2 a[(AVe)*W e ) + {DM* (DM] 



+ cos 2 a[(D*^y(D* k ^) + (DlAT(DM]j ~ 5\C (19) 

Finally, for the generator of the canonical scale transformation we get 



5 C Q -- I d 3 x x k T° k = I d 3 x x k 



E t F kl + ~{P e D k j> e - Dir^e + WDM - D^Vi} 



(20) 
It yields the following gauge covariantized version of ([12]), (|T^), 

5 c s E k = x'diEu + 2E k + x k (diEi + etyWi - ij>*M) (21) 

5 c s A k = x%A k + A k - d^x'A,) (22) 

3 3 

tfj il>e,i = x l diip e4 + -^ ± iex'Aiip^x'di^ + -ip*^ (23) 

In particular, for each of the canonical variable (ip e ,i, i>ti) ^ ne sca ^ e dimension is 3/2 so 
that the canonical scale generator commutes with the number operators for the charged 
particles 

{<S?,JV e)i } = $ f d 3 x i>*e,i<l>e,i = (24) 



We now proceed to inspect the consequences of Shafranov's virial arguments. For 
this we remind that a static minimum energy configuration must be a stationary point of 
the energy ([$]), flI8| ) under any local variation of the fields. Since the scale transformation 
(|P^), ( |13"D is a non-local variation it does not need to leave the energy intact, unless it also 
preserves the pertinent boundary conditions. To determine these boundary conditions, 
we consider the plasma in two different physical environments: 

In the first scenario we have an isolated, localized plasma configuration in an other- 
wise empty space, with a definite number of charged particles 

N e + N t = fd 3 x (^e + VM) (25) 



Since the canonical scale generator commutes with the individual number operators (24) , 
the ensuing variation of the fields is consistent with the boundary condition that the 
number of particles remains intact. By a direct computation we then find for a static 
stationary point of the energy, 



= 5% d 3 x T° n = - / d 3 x T 



i 



r ( 1 E 2 B 2 \ 
J d 3 xi {sin 2 a\D k ^ e \ 2 + cos 2 a\D* k ^ t \ 2 } - ~2 2") ( 26 ) 



i Z7-2 d2 v 

/i 

Since the trace of the spatial stress tensor is a sum of positive definite terms, in analogy 
with Shafranov's virial theorem in ideal magnetohydrodynamics (0) we conclude that 
there can not be any nontrivial stationary points. This means that in an otherwise 



empty space an initially localized plasma configuration can not be confined solely by its 
internal electromagnetic interactions, additional interactions such as gravity must be 
present. Otherwise the canonical scale transformation dilutes the plasma by expanding 
its volume while keeping the number of the charged particles intact, until the collective 
behaviour of the plasma becomes replaced by an individual-particle behaviour of the 
charged constituents. 

The second physical scenario of interest to us describes the bulk properties of a 
plasma: We are interested in an initially localized plasma configuration, located within 
the bulk of a finite density plasma background. In this case the relevant boundary 
condition on the charged fields states, that at large distances their densities approach a 
non- vanishing constant value po which is the density of the uniform background plasma, 

hM 2 ™ /^ (27) 

The canonical scale transformation assigns a non-trivial scale dimension to the charged 
fields. Consequently it can not leave the asymptotic particle density intact, and fails to 
be consistent with the boundary condition fl2"T|) unless po = 0- Instead of the canonical 
version of the scale transformation, we need to employ a non-canonical version of fll2|) . 
(|13|) where the scale dimensions of the charged fields vanish, D^ = 0. When we perform 



the ensuing variation of the fields in the energy density (|i~8|) , instead of (p6|) we find 

~E 2 B 2 1 



d 6 x T u n = d?x 



— — — — (sin 2 a\D k ib e \ 2 - cos^ a I Dtil); 
2 2 2p [ " ! 



(28) 



Now the integrand acquires both positive and negative contributions, which implies that 
a virial argument can not exclude the existence of stable finite energy solitons. Indeed, in 
it has been argued that stable knotted solitons are present. These solitons are formed 
within the bulk of the plasma, in an environment with an asymptotically constant back- 
ground density. A physical example of such an environment is the solar photosphere, the 
solitons are natural candidates for describing stable coronal loops. Another, somewhat 
more hypothetical example could be the ball lightning, in the background of Earth's 
atmosphere. Such solitons could also become relevant in identifying particularly stable 
plasma configurations in fusion experiments, when the plasma is kept at finite density 
by the boundaries of an appropriate vessel. 

We shall now proceed to demonstrate, that the virial theorem ( |2"%| ) is also consistent 
with an appropriate canonical scale tranformation. For this we first notice that ex- 
cluding the kinetic terms, the Lagrangian ( |TTD coincides with that of relativistic scalar 
electrodynamics with two flavors of scalar fields, 

C = \{d, + iA,)<t> x \ 2 + |(^-^)0 2 | 2 - V{4>) - \f 2 u (29) 

Here we have included a Higgs potential V(0), to ensure a non- vanishing asymptotic 
value for the charged fields. For example, we can choose V(<p) oc (</> 2 + §\ — pi) 2 . The 
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Hamiltonian version of (|29|) is 

C = 7r*<9 o 0i + 7Ti<9 o 0t + ^2^)02 + n 2 d (j)* 2 + E^Ai - \(d k + iA k ) 

E 2 B 2 



: -\(d k -iA k )<f)2 



-7T*7Ti — 7T27T2 — V{<j>) 



■A (diEi + Z7T*0l - i-Kl4>\ - 2^202 + ^202) (30) 



2 2 

Notice that now the charged fields are canonically independent variables, a consequence 
of Lorentz invariance. The energy-momentum tensor can be computed directly from 
(|10|). With a proper symmetrization it becomes fully symmetric, as it should since the 
theory is Lorentz invariant. For the energy density we find 



J- r 



I A. 



wi 



+ \D 



kV2\ 



E 2 B 2 
+ ^ + ^ + <*"i + nfa + V((f>) 



Ao{diEi + i(vr*0i - 7Ti0*) - i(n*(p2 - ^fil)} 



The momentum flux is 

T\ = E t F kl + TrJDfc^ + TuDlft + ir*D* k <f> 2 + K 2 D k <p* 2 



(31) 



(32) 



so that instead of ( |23| ) , we find that the canonical scale dimensions of the charged scalar 
fields now vanish. As a consequence the canonical scale transformation is consistent with 
the relevant boundary condition that in the r — > oo limit the system approaches a Higgs 
vacuum <p\ + <p\ — » pjj. This means that the canonical scale transformation must now 
leave the energy of a static stationary point intact which leads to the following virial 
theorem 

= 8% [ d 3 xT° = - I d 3 x r J' : , 



d 3 x 



E 2 B 2 



7r 2 7T 2 ; 



\D k <j>i 



\Dl 



kV2\ 



(33) 



Since the contribution to the pressure from the charged fields is negative, the virial 
theorem can not exclude stable static solitons. 



Finally, we note that even though the static sectors of the two theories (J17Q and (|29f) 
are very similar, these theories actually have a quite different physical content: In the 
relativistic case we may consistently set 7Ti = n 2 = E,i = A = in the static equations 
of motion. This reduces the energy density to a functional form which is manifestly 
magnetohydrodynamical (pi), 



E 



d 3 x 



IZ^I 2 + \D* k fo\ 2 + V{cj>) + 



B 2 



(34) 



But since the canonical scaling dimensions of the charged fields now vanish, the virial 
theorem does not exclude the existence of purely magnetic solitons. On the other hand, 



if in the non-relativistic case (|17j) we set the electric field to vanish, the equations of 
motion become inconsistent unless the electron and ion charge densities are everywhere 
identical. This leads to a contradiction whenever the Hopf invariant is nontrivial ||, 
j|]]. Hence solitons with a nontrivial Hopf invariant are necessarily accompanied with 
a nontrivial electric field. In particular, this means that their properties can not be 
consistently inspected by pure magnetohydrodynamics. 

In conclusion, we have extended Shafranov's virial theorem from ideal magnetohy- 
drodynamics to classical field theory and related it with Derrick's scaling argument. We 
then employed the virial theorem to inspect soliton stability in a realistic field theory 
model of a two component plasma. In line with ideal magnetohydrodynamics, a scaling 
argument reveals that the field theory model does not support stable isolated solitons in 
an otherwise empty space. But the virial theorem does allow for the existence of stable 
solitons within the bulk of the plasma. These solitons are accompanied by a nontriv- 
ial electric field, hence they can not be probed by magnetohydrodynamics alone. We 
suggest that these solitons are relevant in describing coronal loops in solar photosphere, 
maybe even ball lightning in Earth's atmosphere. They might also become useful in the 
design of particularly stable magnetic fusion geometries. 

We thank A. Bondeson, R. Jackiw, S. Nasir, A. Polychronakos and G. Semenoff for 
discussions. 
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